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Abstract 

The problem of spin-dependent transport of electrons through a finite ar- 
ray of quantum dots attached to ID quantum wire (spin gun) for various 
semiconductor materials is studied. The Breit-Fermi term for spin-spin in- 
teraction in the effective Hamiltonian of the device is shown to result in a 
dependence of transmission coefficient on the spin orientation. The differ- 
ence of transmission probabilities for singlet and triplet channels can reach 
few percent for a single quantum dot. For several quantum dots in the array 
due to interference effects it can reach approximately 100% for some energy 
intervals. For the same energy intervals the conductance of the device reaches 
the value ~ 1 in [e 2 /irh] units. As a result a model of the spin-gun which 
transforms the spin-unpolarized electron beam into completely polarized one 
is suggested. 



1 Introduction 

Spintronics is a novel branch of nanoscience which exploits spin properties of elec- 
trons or nuclei instead of charge degrees of freedom pQ. Perhaps the most current 
efforts in designing and manufacturing the spintronics devices are focused now on 
finding novel ways of both generation and utilization of spin-polarized currents. 
This, in particular, includes the study of spin-dependent transport in semiconduc- 
tor nanostructures, which can operate as spin polarizers or spin filters. Although 
spintronics can have a lot of practical applications j2] - [Hj, one of them is most 
interesting and ambitious: the application of electron or nuclear spins to quantum 
information processing and quantum computation jH] - |13j . 

In this paper we suggest a new type of spin polarizer constructed from the finite 
array of semiconducting quantum dots (QD) attached to a quantum wire (QW). 

All dots of the array are identical and each of them carries the spin 1/2. The 
schematic construction of the proposed device is shown for a single dot in Fig.l and 
for the finite array of dots in Fig. 2. 
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Figure 1: Schematic representation of one quantum dot attached to quantum wires. The 
left and right reservoirs are at different chemical potentials fii and fiR. The horizontal 
arrows represent incident, transmitted and reflected electronic wave. The quantum dot 
carries spin s = |. By E s and E t the singlet and triplet states are denoted. By j! the 
electron spin projection S3 = ±^ before and after the scattering by the quantum dot is 
denoted. If the small voltage eV = — is applied a non-equilibrium situation is 
induced and spin-dependent current will flow through the device. 



r 



11 







LEFT RESERVOIR 



RIGHT RESERVOIR 



Figure 2: Schematic representation of spin filter device (spin gun) fabricated from iV 
quantum dots and quantum wires. For details see the caption for Fig.l 



Using the mathematical modeling based on zero-range potentials ^3] we 
reduce the problem of spin-dependent transport in the device to exactly solvable 
scattering problem. In terms of the scattering data, in particular, in terms of spin- 
dependent transmission coefficients through the device, we calculate the most im- 
portant measurable characteristics: the polarization efficiency P{k) and the conduc- 
tance G{k) of the device. 



2 Assumptions, constrains and choice of model 
parameters 

In this section we formulate the physical conditions under which our mathematical 
modeling is relevant. We assume that in the quantum wire connected with a quan- 
tum dots one-mode propagation of electrons in a ballistic regime can be realized. 
For the ballistic propagation of electrons through the QW the de Broglie wavelength 
17J \b = ^2m*k B T nas ^° ^ e mucn greater than both the mean free path l m f p in the 
material of QW and the length d of QW between the neighboring quantum dots, 
i.e. X B > l m fp and X B > d. 

On the other hand in order to use zero-range potentials for mathematical model- 
ing of scattering by quantum dots the de Broglie wavelength has to be greater than 
the mean size ro of the dot: Xb ^ ?"o |14| . There are different ways to manufacture 
quantum dots (see, e.g. [TH], [IH], [20] )• Depending on technology the mean size of 
quantum dots varies from 20 nm for small dots up to more than 100 nm for the large 
dots. In order to satisfy all the constrains mentioned above we have to assume that 
spin-gun can be realized on relatively small QD using narrow-gap semiconductors 
(i.e. semiconductors with small E g — see Table 1). For such materials (see Table 
1) X B is about 50 nm at room temperature and more than 100 nm at the liquid 
nitrogen temperature. 

Table 1.0] 



Semiconductor 


E 9 , eV 
(285 K) 


m* 
m e 


Ag,nm 
(285 K) 


Xb, nm 
(77 K) 


GaAs 


1.430 


0.068 


29.9 


57.7 


InAs 


0.360 


0.022 


52.8 


101.5 


Cd x Hgi- x Te 




x = 0.20 


0.150 


0.013 


68.6 


131.9 


x = 0.30 


0.290 


0.021 


54.0 


103.9 


HgTe 


-0.117 


0.012 


71.4 


137.4 


Zno.15Hgo.s5Te 


0.190 


0.015 


63.9 


122.9 


InSb 




0.014 


66.1 


127.2 



It means that we have either to use in our modeling the parameters d, ro <C 100 nm 
or to assume the working temperature of the device less than T = 77K. Below in 
Section 5 we shall take into account these limitations in our numerical simulations 
of the device operation. 

3 A mathematical model of spin transport through 
a single quantum dot 

In this section we consider a model for scattering of a ID electron on a quantum dot 
attached to the quantum wire. Both ID electron and quantum dot are assigned to 
have spin 1/2 (s e = 1/2, — 1/2 and the possible spin projection on a fixed axis 



are (s 3 ) ej d = ±1/2. We emphasize that in the model in question the quantum dot is 
treated as unstructured uncharged scatterer carrying spin 1/2 as a whole. 

Due to the Breit-Fermi theory the spin-spin interaction in frames of this 
simplest model leads to a singular interaction V s in the model Hamiltonian: V s = 
7 s 5(x) ® (1(1), s(2)}, where x is the distance between the electron and the quan- 
tum dot, 5(x) is the ID delta function, j s is the coupling constant, s(l) and 
s(2) are the electron and the quantum dot spin operators respectively: s(j) = 
(si(j), S2(j), ss(j)), j = 1, 2, Si = |(7j where <Ti, i = 1, 2, 3 are the Pauli matrices and 
(5(1), s(2)) = si(l) <g> si(2) + s 2 (l) ® s 2 (2) + s 3 (l) ® s 3 (2). 

Thus, the Schrodinger equation relevant to the proposed model has the form 



d 2 



dx 2 



h + ls 5(x)®(s(l),s(2)) 



s 3 (l), s 3 (2)) = fc 2 ^(x, s 3 (l), s 3 (2)), (1) 



Here J 4 is the unit operator in the total spin space C 2 (g) C 2 and m* is the effective 
mass of electron in semiconductor material of the quantum wire and ® means the 
tensor product, k 2 = E^-, s(j) = ^ _1 s(j), 7 S = 2m*j s . The wave function \l/ 
depends on the coordinate variable x and s 3 (l), s 3 (2) can take the values ±1/2. 

In accordance with the extensions theory approach jTH] the equation (JTJ) is equiv- 
alent to the following boundary problem 



d 2 



dx 2 



L $ = k 2 ^ 



(2) 



x=+0 



^1 



*'U=-o = 7 a (s(l),s(2))^U=o- 



(3) 



In order to solve the problem (j2J, (jSJ) one has to separate the spin variables, i.e. 
to calculate the matrix elements of the Hamiltonian from equation ([IJ in the spin 
basis constructed below. Let us denote by x±(l) an d X±(2) the spin part of the 
electron and the dot wave functions respectively with the value of spin projection 

s s(j) = = tl/2 ) J — 1)2. One can take these spinors in the form x+(J) ~ 







X-(j) 




It is well known |23 j that the two-spin space Ti Sl ® H S2 ~ C 2 ® C 2 



can be decomposed into the direct sum of two invariant (with respect to operators 
of representation of the group of rotation) subspaces: H Sl <%> TC S2 — D © Dx, where 
-Do is the singlet space (subspace of the total spin 5 = 0) which is spanned on the 
spinor Xo,o, dimD = 1, and D% is the triplet space (subspace of the total spin 
5=1) which is spanned on three spinors {xi,i> Xi,o, dimDi = 3. 

The orthonormal basis of spinors in the triplet subspace has the form j2H]: = 
X+(l)®X+(2), = X-(l)®X-(2), Xi.o = ^(x+(l)®X-(2)±X-(l)®X+(2)). The 
second index in Xi,s 3 means the values of the total spin projection 5 3 = 1, 0, —1 The 
basis vector in the singlet subspace reads Xo,o — 7^(x+(l) ®X-(2) — X-(l) ( 8 ) X+(2))- 
The above described spinors form the orthonormal basis in the total spin space 



TC S1 <8> H S2 ~ C 2 <8> C 2 . Hence the total wave function s 3 (l), s 3 (2)) can be 
decomposed as follows 



*(x, s 3 (l), s 3 (2)) = ^i,i(x)xi,i + ^i,o(^)Xi,o + + ^o,o(»Xo,o (4) 

In order to separate the spin variables in the problem (|5J) it is sufficient to 
calculate the action of the operator (s(l),s(2)) on these spinors. One can easily 
verify that (s(l), s(2)) xi,s 3 = zXi,s 3 , where 5 3 = 1,0,-1 and (s(l), s(2)) Xo,o = 
— 4X0,0- 

Since the chosen basis of spinors is an orthonormal system the separation of 
spin variables in the channel 5 = and 5 = 1 gives the following boundary value 
problem for the coordinate parts of the wave function (@J) 

<h I 2 - k 2 ) if) =0 (5) 



dx 2 

if)(+0) = v<-o) 
fUo-^=-o=( 7 U 7 ° oo ) V(0). (6) 

Here if) is a two-component vector if) = \ }~' S3 ] (x) where 5 3 = 1, 0, —1, 700 = — f 7s 

is the coupling constant in the singlet channel and 711 = is the coupling constant 
in the triplet channel. In the proposed model in the absence of relativistic effects 
the total spin 5 conserves and the triplet-singlet and singlet-triplet transitions are 
forbidden. As a consequence, the antidiagonal elements of the coupling matrix in 
boundary condition © are equal to zero. Since the projection of the total spin 5 3 
also conserves, one can choose in the boundary problem above any value of this 
projection. In what follows we set 5 3 = 1. 



4 Scattering on periodic array of N quantum dots 

The single-dot Hamiltonian (JJJ) for the case of periodic array of iV quantum dots is 
modified as follows 

h 2 J 2 l=N 
H {N) = -^-^ ® h + E ls5{x - Vl ) ® 5(2)}, (7) 

where points yi are separated by distance d. 

All the constructions of the previous section are easily generalized for this case. 
Hence after separation of spin variables in the equation if( 7V )v|/( Ar ) = _E\l>( Ar ) one 

obtains for the coordinate part if)^ N '(x) = ( ,(]v) I i x ) °f the total wave function 



vpW the following boundary value problem 



'0, 



0,0 



^ 2 ®h-k 2 )^ =0 (8) 



V> (JV) (w + 0) = ^W(y,-0) 



In 
7oo 



^ N \ Vl ), l = l,2,...,N. 



(9) 



The transmission coefficient in the triplet channel T^\k) and singlet channel 
T o W are calculated in terms of elements of the 4x4 matrix T = j = 

1,2,3,4, (i is the line number, j is the column number) which is the product of 
transition matrices through dots and wires connected dots. It reads T = Q N ~ 1 Qi, 
where 



Q 2ik 



( 7n + 2%k 

7oo + 2ik 



e 2tdk j n 







V 



-7ii 






-7oo 



-2idk 





{2ik- ln ) 




e 2idfc 7oo 



\ 





(2ik - 700) 7 



(10) 



The matrix Qi differs from Q by absence of multiplier e~ 2ldk in its matrix elements. 
Here d is the distance between the centers y\ where 5— functions are localized. As 
applied to the physical situation we are modeling d is the length of a quantum 
wire connecting two neighbor dots. The transmission coefficients have the form 
Tip = t n + t 13 e- 2ikyi R[i \ = £ 22 + h^^R^. The reflection coefficients 



RyP in the triplet and singlet channels respectively are given by equations 
= e 2lk yW-\t 41 t u - £31*44), R$ ] = e 2lk yW-\t 43 t 32 - £33*42), where D = 
*33£44 — £34*43. The relations \T^\ 2 + \R^\ 2 = 1 (j = 0,1) which imply the 
unitarity of scattering has been checked by means of analytical calculations package. 



In the case N 
(j = 0, 1) 



1 one has = 2ik(2ik — jjj) 



(i) 



1 R ( 



5 Results and discussions 

In this section we present the results of the numerical calculations of the operation 
regime characteristics for the spin gun and discuss the proper choice of the model 
parameters as well as the choice of semiconductor materials for QW and QD. We 
use the applied bias V : eV = hl — Hr for tuning the spin polarization effects in 
the device. 

In order to calculate the conductance G of the device we exploit the two-channel 
Landauer formula at zero-temperature limit: 



e 2 

G = h 

irh 



T$\E F + eV)\ 2 + \TlP(E F + eV)f 



where Ep is the Fermi level of QW. 

For unpolarized incident beam of electrons the scattering by QD array in the 
triplet and singlet channels leads to different transmission coefficients and Tqq\ 
Following we shall use the polarization efficiency 



which determines the difference of transmission probabilities through the spin gun 
for the triplet and singlet spin states. 

At a given wave number k the conductance G and polarization efficiency P of 
the device both depend on the model parameters: the distance d, the coupling 
constants 711, 700 and the number N of QD. They also both depend on the choice 
of semiconductor materials from which the QD and QW are fabricated, i.e. effective 
mass to*, Fermi energy Ep and mean size r of quantum dots. Varying the bias 
applied one can calculate the dependence of G and P on k = \J 2m (Ip +IZ) a f fixed 
model parameters and chosen to*, Ep and r . Let us mention that from spin analysis 
in Section 3 one has 700 = — 37n. Since d and N are free parameters of the model 
the only parameter has to be fixed is 711. The reasonable way to fix this parameter 
is to take it proportional to the inverse mean size of the quantum dot. For relatively 
small QD r ~ 20 -j- 30 nm. Hence 711 ~ 0.03 -j- 0.05nTO _1 . Taking d ~ 5 Own, we 
obtain the value of order £s 1 for the dimensionless parameter 7 n = 711 -. 

In Figures 3 - 5 we show the dependence of G and P on the dimensionless 
parameter k = k^ for 7 n | = 1 in the cases N = 1,2,8. We also indicate the 
position of Fermi level E F = 35meV in InSb (to* = 0.014to 6 ) and fix the inter-dot 
distance d — 45 nm. 

The general behavior of the conductance G and polarization efficiency P for the 
wide range of parameters variation is the following. On the fc-axis with the growth 
of the number N of quantum dots there arise the so-called "working windows" (see 
Figure 5), i.e. such intervals [k\, /C2], [^3, k^],...m which P and G are sufficiently 
high simultaneously. Due to the interference processes in a QD array at N ~ 8 -j- 10 

2 

the values of P and G on the working windows approach P 100% and G ~ \ 
which equals to one half of its upper limit. It means that at k lying in the working 
windows the spin gun really produces the fully polarized spin beam and conducts 
at fairly high level. In Figures 3-5 the working windows lie on the right of Ep up 
to the value of k — ^\J 2m * <yE J ^ +eV ^ which corresponds to the bias applied. 

6 Plans and prospect 

There are a lot of interesting modifications of the model we constructed. For ex- 
ample, we are going to take into account the internal structure of quantum dots. 
We plan to supply the quantum dot by prescribed internal structure remaining the 
model exactly solvable. Other interesting modification is to take into account the 
Coulomb interaction between the electron and quantum dot if the dot is charged. 
We shall prove that if we switch on the Coulomb interaction the modified model 
remains exactly solvable. Finally we are going to incorporate in our model the 
magnetic field in order manipulate spins in QW and QD. 
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kdfn [dimensionless] 

Figure 3: The conductance G and polarization efficiency P versus dimensionless momen- 
tum k = k^ of incident electrons in the case of one quantum dot. By the dashed line the 
Fermi level Ep = 35 meV in InSb is indicated. The interdot distance d = 45 nm. 




kdM [dimensionless] 



N =2, y 11 d/7t =1 , 7 00 d/ir =-3 




kdM [dimensionless] 

Figure 4: The conductance G and polarization efficiency P versus dimensionless momen- 
tum k = k^ of incident electrons in the case of two quantum dots. By the dashed line the 
Fermi level Ep = 35 meV in InSb is indicated. The interdot distance d = 45 nm. 




kdM [dimensionless] 



N=8, T„d/7i=1, Y 00 d/ic=-3 




taMi [dimensionless] 

Figure 5: The conductance G and polarization efficiency P versus dimensionless momen- 
tum k = k^ of incident electrons in the case of 8 quantum dots. By the dashed line the 
Fermi level Ep = 35 meV in InSb is indicated. The interdot distance d = 45 ran. 



